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Abstract. L1521E seems unique among starless cores. It stands out in a distribution of a ratio (R) that we define to asses core 
evolution, and which compares the emission of the easily-depleted C 18 molecule with that of the hard to deplete, late- time 
species N 2 H + . While all cores we have studied so far have R ratio lower than 1, L1521E has an R value of 3.4, which is 8 
times the mean of the other cores. To understand this difference, we have modeled the C' s O and N2H + abundance profiles in 
L1521E using a density distribution derived from 1.2mm continuum data. Our model shows that the C ls O emission in this core 
is consistent with constant abundance, and this makes L1521E the first core with no C ls O depletion. Our model also derives 
an unusually low N 2 H + abundance. These two chemical peculiarities suggest that L1521E has contracted to its present density 
very recently, and it is therefore an extremely young starless core. Comparing our derived abundances with a chemical model, 
we estimate a tentative age of < 1.5 x 10 5 yr, which is too short for ambipolar diffusion models. 
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1. Introduction 

Observations of starless dense cores over the last few years 
have shown that molecular depletion onto cold dust grains 
is a common featur e in low-mass star-fo rming regions (e.g., 
ICaselli et all 1 19991 iTafaUa et all 120021) . As a core con- 
tracts, the amount of depletion of certain molecules like CS 
and CO increases with time ("e.g.. iBergfn & Langerl 1 19971 
lAikawa et al. Il2003l> . so we can expect that depletion will pro- 
vide an accurate clock for core evolution when its details are 
well understood. Even before this is a reality, depletion can 
be used as a qualitative time marker, in the sense that evolved 
cores should be more depleted of certain species than younger 
cores. When we attempt to build a sequence of cores having dif- 
ferent amounts of depletion, however, we find that the popula- 
tion of undepleted (or lightly) depleted cores is missing. These 
undepleted cores should represent the earliest stages of core 
contraction, and probably they have not been identified so far 
because of an observational bias. Given the interest of young 
cores for studies of core contraction, finding a member of this 
population is an urgent challenge. 

To identify cores with low degree of depletion, we have 
started a systemati c study of cores with weak NH3 lines in the 
survey of iBenson & Mversl Jl989l) . as previous observations 
of cores bright in NH3 h ave only found cases of s trong de- 
pletion (e.g.. lCasellietal.lll999tlTafallaetal.ll2002l) . Among 
our targets, the L1521E core in Taurus has revealed itself as a 
core with negligible depletion, and here we report a prelimi- 
nary analysis of its properties. This core has alrea dy been iden- 
tified as a very young core bv lHirota et all j2002), who found 
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it very prominen t in carbon-chain m olecules, comparable to 
TMC-1 (see also lAikawa et~Z1 Eo03). In this paper we show 
that L1521E has the lowest measured level of C ls O depletion, 
and in a forthcoming article we will present the results of an 
ongoing chemical survey of this core. 

2. Observations 

We observed L1521E in C 18 O(l-0) and N 2 H + (l-0) with the 
FCRAO 14m telescope during 2002 March. The observations 
were done with the SEQUOIA array in on-the-fly position 
switching mode. The facility correlator provided a velocity res- 
olution of about 0.07 km s _1 , and all data were converted into 
the main beam brightness scale assuming an efficiency of 0.55. 
The telescope beam size was approximately 50". 

We observed L1521E with the IRAM 30m telescope in 
2003 September (lines) and October (continuum). The line ob- 
servations consisted of maps in N2H + (l-0) and C ls O(2-l), and 
crosses in C I7 O(l-0), C I7 0(2-1), and C 34 S(2-1). All data were 
taken in frequency switching mode with a velocity resolution 
of 0.03 km s _1 , and were converted into main beam brightness 
units using standard efficiencies. The 1 .2mm continuum obser- 
vations were done with MAMBO-2 in on the fly mode. Sky 
dips were used to measure the atmospheric extinction, and ob- 
servations of planets were used for calibration. The resolution 
of the 30m telescope varies from 26" for N2H + (l-0) to 11" for 
C ls O(2-l) and the 1.2mm continuum. 

3. The R ratio and the uniqueness of L1 521 E 

To search for cores at the earliest stages of contraction, it is 
convenient to define a simple observable that measures deple- 
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Fig. 1. Left: R parameter (C 18 O(l-0)/N 2 H + (l-0) ratio at the core center) as a function of the central N 2 H + (l-0) integrated 
intensity. Each point represents a core and the dashed line at R = 1 represents the average R value for the core models with 
no C 18 depletion in lTafalla et all J2002I) . The squares represent cores studied in that paper (L1400K, L1517B, L1498, L1495, 
and L1544 by increasing N 2 H + intensity), the triangles represent cores from an unpublished survey (L981-1, LI 197, LI 155C1S, 
L1155C1N, L1512, L234E-1, L1333, and L694-2, same order), and the star symbol indicates B335, the only core with a YSO. 
Note how L1521E stands out from the core sample. Right: Large scale C 18 O(l-0) integrated intensity map overlayed on a DSS 
blue image. Contours are .2, .4, ... K km s _1 . All radio data are from FCRAO. 



tion without the need of a complicated radiative transfer cal- 
culation. One possibility is to compare the emission of two 
molecules known to have different depletion behaviors at typ- 
ical core densities, and for this we choose C 18 and N 2 H + . 
C 18 is a heavy depleter at densities of a few 10 4 cm 4 , while 
N 2 H + seems to remain in the gas phase up to densities o f about 
10 6 cm" 3 (e.g.. lCaseUi et al. ll999ilTkfalla et al. lEooll) . Thus, 
if we define R [C ls O/N 2 H + ] (R hereafter) as the ratio between 
the integrated intensities of C 18 and N 2 H + at the center of a 
core, we expect R to be a sensitive indicator of depletion. In 
particular, we expect R to be low in cores with strong C 18 
depletion and high in cores with little C 18 depletion. 

Our definition of R has the additional advantage that it com- 
pares two molecules with different gas-phase formation time 
scales: C ls O forms relatively quic kly in the core histor y while 
N 2 H + is a late-time molecule (e.g.. lAikawa et al . 2003). Thus, 
if chemical youth is correlated with lack of depletion, the R ra- 
tio should be further enhanced in young cores because of their 
relatively lower N 2 H + abundances. A signature of youth and 
lack of depletion, therefore, is the combination of a large R and 
a relatively low N 2 H + intensity. 

To calibrate the R ratio, we test its behavior in a sample 
of cores selected for their strong N 2 H + or NH3 emission, and 
therefore likel y to be evolved . This sample includes the 5 cores 
studied in Tafal la et al. I faOOA a set of 8 additional cores from 
an unpublished survey, and a core with an embedded YSO, 
B335. All these cores have been mapped in C ls O and N 2 H + 
with the FCRAO telescope, and their R ratio has been esti- 
mated at the position of maximum N 2 H + emission. In addi- 
tion and as a reference, we have estimated the R ratio for the 
5 Monte Carlo core models with constant C 18 abundance in 
Tafal la et al. I J2002h . finding an average R value of 1 .0. 

As Figure 1 illustrates, all cores in our sample lie below the 
R — 1 line despite their factor of 4 spread in N 2 H + (l-0) inte- 
grated intensity. The sample has a mean R value of 0.43 with 
an rms of 0.14, and its location under the R = 1 line is sug- 
gestive of strong C ls O depletion (a conclusion supported by 
detailed radiative transfer modeling of some objects). This sys- 
tematic behavior suggests that R = 1 can be taken as a reason- 



able boundary between cores with and without C O depletion, 
and that a search for undepleted cores should look for objects 
above the R = 1 line. Thus, over the last few years we have car- 
ried out a search for high-/? cores using the FCRAO telescope, 
mapping core candidates both in N 2 H + (l-0) and C 18 O(l-0). 
Although the search is still ongoing and several candidates 
seem promising, one object stands out because it combines 
the highest R value and the lowest N 2 H + intensity. This core, 
L1521E, shown in Fig. 1 as the outlying point labeled with its 
name, has R = 3.4, which is 8 times larger than the mean or 
20 times the rms. These characteristics make L1521E the best 
candidate we know of to an undepleted core, and as we will in 
the following section, a detailed analysis of its molecular com- 
position confirms this interpretation. 

4. The internal structure of L1521E 

The high R ratio of L1521E can arise from a high central C ls O 
abundance, a low N 2 H + fraction, a combination of both, or 
some excitation effect (after all, the two molecules have very 
different dipole moments). To discern between these options, 
we now solve the radiative transfer of the two species and deter- 
mine their abundance profiles. For this, we rely on observations 
of the 1.2mm continuum, C 18 O(l-0), C ls O(2-l), C 17 O(l-0), 
C 17 0(2-l), and N 2 H + (l-0), shown in part in Figure 2. 

4.1. Physical parameters 

We first derive the density profile of LI 52 IE from its distribu- 
tion of the 1.2mm continuum emission. This emission is op- 
tically thin and appears unaffected by chemical differentiation 
or excitation effects, so it can be easily inverted into a den- 
sity distribution (e.g.. lAndre et al. Il996l) . As Fig. 2 shows, the 
1.2mm emission (like that of the molecules) is elongated NW- 
SE, and therefore deviates from spherical symmetry. Thus, to 
derive a density distribution that can be used by our spherically- 
symmetric Monte Carlo radiative transfer analysis (see below), 
we first average elliptically the observed emission, assuming 
an aspect ratio of 2:1 and a position angle of -60°. Then, 
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Fig. 2. Maps of L1521E in 1.2mm continuum, C 18 O(l-0), C 18 0(2-l), N 2 H + (l-0), and C 34 S(2-1). Note the similar distribution 
of all tracers, something not seen in cores with depletion. First contour and contour spacing are 4 mJy beam -1 for the 1.2mm 
continuum map, 0.2, 0.3, 0.1, and 0.1 K km s" 1 for C 18 O(l-0), C ls O(2-l), N 2 H + (l-0), and C 34 S(2-1), respectively. Offsets are 
referred to a 2 ooo = 4 h 29 m 15 s .7, fcooo = 26°14'5". Note that the C 34 S(2-1) map is heavily undersampled. The C 18 O(l-0) map is 
from FCRAO data and the other maps are from IRAM 30m data. 



we find an analytic density profile that after simulating an on- 
the-fly observation predic ts the correct radial p rofile of 1 .2mm 
continuum emission (see Tafall a et al. 1120021 for details). For 
consistency with previous work, we assume a dust emissiv- 
ity K\2mm = 0.005 cm 2 g -1 and a dust temperature of 10 K 
(see below). In this way, we derive a density distribution of 
the form «(r) = «o/(l + ( r / r o) a )> where «o - 2.7 x 10 5 cm -3 , 
r = 6.3 x 10 16 cm (= 30" at 140 pc), and a = 2. Fig. 3 shows 
that this profile provides a reasonable fit to the observed 1 .2mm 
continuum emission. 

The NH^ emission from L1521E is extremely weak 
feenson&Mverslll989t iHirotaet al.ll2002h . so we rely on 
C I7 to determine the gas kinetic temperature. C 17 seems un- 
depleted at the core center (see below), and fits to its hyperfine 
structure indicate that the emission is optically thin. Thus, we 
can convert the observed central 2-1/1-0 intensity ratio (« 1.3) 
into an excitation temperature (« 9 K) which we approximate 
as 10 K. The 2-1/1-0 line ratio seems constant with radius, 
so the excitation temperature must also be close to constant. 
Given the low dipole moment of C 17 (0.1 D), this excitation 
temperature should be a good approximation to the gas kinetic 
temperature, and we will assume in our modeling a constant 
gas temperature of 10 K. 

Finally, we characterize the velocity field of L1521E using 
hyperfine fits to its N2H + (l-0) emission. These fits indicate an 
intrinsic N 2 H + linewidth of 0.3 km s _1 , which means a nonther- 
mal component of about 0.27 km s _I (FWHM). When model- 
ing the C 18 O(l-0) lines, a slightly larger linewidth is needed 
in the outer core, and we have parametrized this behavior as a 
linear increase in the non thermal component up to 0.58 km s _1 
at 4 x 10 17 cm (the core outer radius in our model). 

4.2. Molecular abundance profiles 

We solve the radiative transfer for N 2 H + and the CO iso- 
topomers to derive their rad ial profiles o f abunda nce. For this, 
we use a Monte Carlo code (Bernes 1979; Tafalla et al. 120021 
120041) and fit simultaneously both the radial profiles of inte- 
grated intensity and the emergent spectra of all transitions (the 
radial profiles have been derived with the same elliptical aver- 
age as used for the continuum). Given the previous parameter- 



ization of the LI 52 IE core, the only free parameters to fit the 
data are the abundances of each species. 

To determine the N 2 H + abundance we use the 1-0 tran- 
sition. This transition has 7 hyperfine components, and our 
model treats each component as a separate line, taking explic- 
itly into a ccount the hyp erfine splitting of the first 3 rotational 
levels (see Tafalla et al. 2004 for further details). We have run 
this model assuming constant N 2 H + abundance, and as Figure 
3 shows, the predicted emission for a value of 2 x 10~" fits 
simultaneously the radial profile of integrated intensity and 
the central emerging spectrum. A constant N 2 H + abundance is 
also found in other cores, but the value we derive for LI 52 IE 
is about 8 times lower than that d erived for the L1498 and 
L1517B cores using the same code JTafalla et al. 1120 04). This 
difference is highly significant and makes L1521E the star- 
less core with the lowest N?H + abundance ( note that the N 2 H + 
abundances estimated by Ta falla et al. l2002l are artificially low 
because of their simplified treatment of the hyperfine structure). 
The low N 2 H + abundance explains the location of L1521E in 
the far left of Fig. 1 . 

To determine the abundance of the CO isotopomers, we 
now fit the radial profiles and emergent spectra of C 18 0(1- 
0), C I8 0(2-1), C I7 O(l-0), and C n O(2-l). We firs t fix the 
C 18 Q/ C 17 abundance ratio to its ISM value (3.65. IPenziasI 
1981), so we are left with one abundance profile as free param- 
eter. In contrast with what is found in other cores, the LI 52 IE 
data can be fit with a model of constant CO abundance. This 
is illustrated in Fig. 3 and makes L1521E the first core with no 
sign of CO depletion. Our best-fit abundance is X(C 18 0) = 0.7 
xl0~ 7 , very similar to that found in the outer parts of L1498 
and L1517B, althoug h about 2.5 times low er than the stan- 
dard determination bv lFrerkingetal.Nl982l) . This lower value 
may indicate a slight error in the dust parameters or that the 
iFrerking et al. I dl982h abundances need revision. In any case, 
it is the combination of high CO and low N 2 H + abundances 
what causes the hight R ratio of L 1 52 1 E. 

Although an analysis of the CS emission in L1521E is out 
of the scope of this paper, we note that preliminary models indi- 
cate that this species also has constant abundance all the way to 
the core center (already suggested by the C 34 S(2-1) emission 
in Fig. 2). CO and CS, therefore, seem not depleted in L1521E. 
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Fig. 3. Radial profiles of integrated intensity (left) and central emerging spectra (right) for 1.2mm continuum, N2H + (l-0), 
C 18 O(l-0), C ls O(2-l), C 17 O(l-0), and C I7 0(2-1) (data origin as in Fig. 2). The squares and histograms are observed data, 
and the solid lines are model predictions. The model fitting the 1.2mm continuum is the basis of our core density profile, and the 
models fitting the line emission have been calculated using a Monte Carlo radiative transfer model. All models assume constant 
molecular abundance (see text). 



5. An extremely young core 

The absence of CO (and CS) depletion in LI 52 IE seems not to 
result from molecules finding it harder to freeze out in this core. 
Its central density is typical, as well as is its gas temperature of 
10 K (the dust temperature cannot be much higher given the 
observed central density). The very low abundance of the late- 
time molecule N2LF further suggests that the lack of depletion 
in LI 52 IE results from this core being less ch emically pro- 
cessed than others, as it h a s also been argued by Hi rota et al. I 
J2002l) and lAikawa et al. I d2003l) . L1521E, in fact, combines 
the lowest C I8 depletion and lowest N2H + abundance, so it 
appears to be the less processed and therefore youngest core 
known. 

To estimate how young LI 52 IE is, we compare our ob- 
servations wi t h the predictions from the recent model of 
Aika wa et al. ] J2003), who have studied the evolution of a 
core undergoing Larson-Penston contraction taking into ac- 
count both gas-phase and grain-surface chemical reactions. 
Although the exact form of contraction is unclear, this model 
seems preferable because it covers the chemical evolution of 
the gas since time prior to molecular depletion when it has an 
initial density of 10 4 cm -3 . Thus, we have taken the model 
abundance profiles at the earliest time calculated by the au- 
thors (1.52 x 10 5 yr, their Fig. 2a) and predicted intensities 
for the conditions of L1521E. In addition of having to multi- 
ply the predicted abundances by a; 3 to approximately fit the 
data, we find that for C 18 the slight central drop predicted by 
the chemical model is barely consistent with the data, and that 
the predicted central increase of N2FF abundance is already 
too pronounced (the central density at this early time is close 
to our estimate). This indicates that the model predicts more 
chemically processing at time 1 .52 x 10 5 yr than observed, even 
for the relatively fast Larson-Penston flow. With the caveat that 
the time estimate depends sensitively on the assumed sticking 
probabilities, we conclude that the comparison implies a con- 
traction age of < 1.5 x 10 5 yr. 



Our contraction age for LI 52 IE seems consistent with the 
estimate bv iLee & Mvers (1999) that the typical lifetime of 
starless cores is 0.3-1.6 Myr, as this core is younger than oth- 
ers. This age, however, seems incompatible with models of core 
formation via ambipolar diffusion, as they require evolutionary 
times of the order of 1 My r in the most favorable conditions 
fe.g.. lCiorek&Basul200ll) . Thus, if L1521Ehas significantly 
contracted on a time scale as short as suggested by chemical 
models, a more dynamical picture seems needed. The colli- 
sion of supersonic flo ws, as proposed by the supe rsonic tur- 
bulence scenario (e.g.. iMac Low & KlessenlEo 03). however, 
seems also inconsistent with the observations, given that the 
turbulent component of the linewidth {cr~0.ll km s~') is still 
subsonic, and that the velocity differences across the L1521E 
core are only of the order of the sound speed (a description 
of the core velocity structure will presented in a future paper). 
Clearly further analysis of LI 52 IE is needed to understand its 
apparent rapid formation and lack of molecular depletion. At 
this point, it is suggestive to speculate that LI 52 IE belongs to 
an up to now overlooked population of extremely young cores. 
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